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ON A SOLUTION OF THE NONLINEAR DIFFERENTIAL EQUATION FOR TRANSONIC FLOW 

PAST A WAVE-SHAPED WALL 1 

By Carl Kaplan 


SUMMARY 

The Prandtl-Busemann small-perturbation method is utilized 
to obtain the flow of a compressible fluid past an infinitely 
long wave-shaped wall. When the essential assumption for 
transonic flow ( that all Mach numbers in the region of flow are 
nearly unity) is introduced , the expression for the velocity poten- 
tial takes the form of a power series in the transonic similarity 
parameter. On the basis of this form of the solution, an attempt 
is made to solve the nonlinear differential equation for transonic 
flow past the wavy wall. The analysis utilized exhibits clearly 
the difficulties inherent in nonlinear-flow problems. The inves- 
tigation has nevertheless been rigorously carried to the point 
where the question of the existence or nonexistence of a mixed 
potential flow free of discontinuities can be settled by the behavior 
of a single power series in the transonic similarity parameter. 
The calculation of the coefficients of this dominant power series 
has been reduced to a routine computing problem by means of 
recursion formulas resulting from the solution of the differential 
equation and the boundary condition at the wall. One of the 
interesting results of the analysis is the rigorous statement that 
the transonic similarity parameter must be less than four-thirds. 

INTRODUCTION 

The present report considers in detail a notion first ex- 
pressed in reference 1, namely, that the Von Kfirm&n transonic 
similarity rule is implicitly contained in the potential-flow 
solution for high-subsonic flow past a prescribed body. The 
calculations are very much simplified by the choice of rec- 
tangular Cartesian coordinates as independent variables. For 
this reason the solid boundary chosen is a two-dimensional 
wavy wall of small amplitude extending to infinity in 
both the downstream and upstream directions. The problem 
is first treated by means of the Prandtl-Busemann iteration 
method for high-subsonic undisturbed speeds. On the basis of 
tins solution an attempt is made to solve the nonlinear differen- 
tial equation for transonic flow (corresponding to the linear 
Tricomi equation in the hodograph plane) past a wavy wall. 
The analysis utilized exhibits clearly the difficulties inherent in 
nonlinear-flow problems. The investigation has nevertheless 
been rigorously carried to the point where the question of the 
existence or nonexistence of a mixed potential flow free of 
discontinuities can be settled by the behavior of a single power 
series in the transonic similarity parameter. The calculation 


of the coefficients of this dominant power series has been 
reduced to a routine computing problem by means of recur- 
sion formulas resulting from the solution of the differential 
equation and the boundary condition at the wavy wall. 

The author wishes to acknowledge the invaluable aid and 
advice of Dr. A. Busemann of the Langley Laboratory during 
the writing of this report, especially with regard to the final 
section “General Analysis.” 

CALCULATION OF HIGH-SUBSONIC FLOW PAST A WAVY 
WALL BY MEANS OF PRANDTL-BUSEMANN 
ITERATION METHOD 

The fundamental nonlinear differential equation for the 
potential flow of a compressible fluid can be written as 

“ 2 2 M 4>y Y 2 M affxfyfxY = 0 

( 1 ) 

where 

£- 2 =l+*=^ MJ[1-{u 2 +v 2 )\ (2) 

and 

f velocity potential of flow 

X, Y rectangular Cartesian coordinates in plane of flow 
u, v fluid velocity components along X- and F-axis, 
respectively 

U undisturbed stream velocity 

c local speed of sound 

speed of sound in undisturbed fluid 

M„ Mach number of undisturbed stream (U/cf) 

y ratio of specific heats at constant pressure and con- 

stant volume 

The quantities f, X, Y, u, and v are nondimensional with a 
characteristic length l as unit of length and the undisturbed 
stream velocity TJ as unit of velocity. The subscripts X and 
F in equation (1) denote partial differentiation with respect 
to the designated variables. 

In order to obtain the Prandtl-Busemann iteration equa- 
tions based on small perturbations of the undisturbed 
stream, the assumption is made that the velocity potential 
<}> can be expanded in the form 


<t> — X-fftYfiYfsY ... (3) 

1 Supersedes NACA TN 2383, “On a Solution of the Nonlinear Differential Equation for Transonic Flow Past a Wave-Shaped Wall” by Carl Kaplan, 1951. 
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For the purpose of defining and controlling the iteration 
procedure, the function <j> n +i and its derivatives are re- 
garded as small compared with the preceding approxima- 
tion <p n and its derivatives. From equation (3) and the 
fact that for irrotational or potential flow u—<j) X and v=<t>r, 


and 


'W=l+0IX+<f>2A' + <f>3X + • • • 

V=4>iy + <f> 2 r + <f>3r + • ■ • 


When these expressions for u and v are introduced into 
equation (1), together with the corresponding expression 
for c 2 /cj given by equation (2), and the powers and prod- 
ucts of < fi n and their derivatives are grouped according to the 
assumptions of the small-perturbation method, the following 
iteration equations for the first three approximations result : 

$lxx + 0 lj/y : — 0 (4) 

4 ) 2xx~\~ 4 > 2yy^ = 2 [(1 + &')<{)lx4 > lzz~\- ^ly^lxyl (o) 

<j>3xz J r4 l 3v«—2 M „ 2 j(l + cr) (4’lxx4>2x~ l r 4>lor<l>2i:x)-\- 
[2/3 2 (l + <r) — 1] 4>lx<f>ly^lxy J r 

(1 4- 0-)|jf/3 2 (l + a) ~^4>ixx<t>lx i +-^(o-p 2 — l)4>lxx4>ly 2 + 

< Plxi/ < p2y~\~ <t>ly4 > 2xy ( ( 6 ) 


where x and y are new independent variables defined by the 
transformation 


and 


y=pY 


/3 2 =1 ~MJ 

__Y + 1 MJ 
a 2 B 2 


(7) 


Equation (4) is a Laplace equation and equations (5) and 
(6) are Poisson equations where the right-hand sides contain 
only previously determined quantities. These equations 
have been treated recently in reference 2 where the particu- 
lar integrals of equations (5) and (6) are given in real form. 
These particular integrals have been utilized in obtaining 
the flow over the wavy wall. Thus, the equation of the 
infinitely long wave-shaped wall (fig. 1) is assumed to be 

Y=a cos aX (8) 

where 

a amplitude of wave 



X wave length 


r 



Figure 1. — Wave-shaped wall. 


The reference element of length is conveniently chosen to 



The equation of the wavy wall in nondimen- 


sional form and in terms of the variables x and y then becomes 


y= j3e cos x 


(9) 


where, if the thickness coefficient of the wavy wall is defined 
as > then e=f. 

A/ ITT 

The expression for the velocity potential 4>, obtained with 
the aid of the particular integrals given in reference 2, 
including terms of the order e 3 and satisfying the boundary 
conditions at infinity and at the solid wall to the same order 
is as follows: 


<t>=x+^e v sin { 2 <rM to 2 2 / + [2 + <r + 

(2 — o’)/? 2 ] } e~- v sin 2x—~ ({ -8 + M ro 2 [24 + 2«r- 

llcr 2 — cr(22+13cr)/3 2 ]} e~ v -\-M„ 2 [<x (6 + 5ff) + 

(8 + 10 a + 3 a 2 ) P 2 \e~ 3u + 4: a (2 + a) MYye~ 3v ) sin x + 

~(£f {{27 + MJ (- 25 + 11 <7 + 4 < 7 2 -9/3 2 )]e- 3 ‘'+ 

QMJa(3 + 2c + Sp 2 )ye- iv +9a 2 M a Ye~ 3v ] sin 3z+ . . . 

( 10 ) 

It is of interest to examine equation (10) when the as- 
sumptions of transonic-flow theory are introduced. These 
assumptions are essentially that the undisturbed flow veloc- 
ity differs only slightly from the speed of sound, that the 
velocity component normal to the oncoming flow is small 
compared with the speed of sound, and that the velocity 
component in the direction of the oncoming flow is of the 
order of the critical velocity c*. If the undisturbed stream 
is in the direction of the positive x-axis, then the velocity 
potential 4>, referred to the critical velocity c*, can be written 
as (see reference 1) 

<f>=a: + ^Y (1 — MJ)j{x,y) (11) 
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The second term on the right-hand side of this equation 
is the disturbance-velocity potential and implies that terms 
involving powers of 1 —MJ higher than the first are to be 
neglected. The differential equation satisfied by the func- 
tion f(x,y) is obtained from equation (1) and takes the fol- 
lowing simplified nonlinear form: 


by 2 bx bx 2 


( 12 ) 


The boundary conditions to be fulfilled by f(x,y ) are as 
follows: 


If terms of only the first power of l—MJ are retained, it 
follows from equation (15) that 

M*= 1 ]-r (1-MJ+ . . . (16) 

y-r i 

Then, multiplying the right-hand side of equation (10) by 

this expression for M*, replacing e//3 by — (1 —M„ 2 )k 

from equation (14), and neglecting all terms containing 
1 — M„ 2 to a power higher than the first yields the following 
expression for <j>: 


1 1 



bx 



J 

► 

(at y= co) 

bf , . 

x~— — k sm x 
by 

(at y=0, - 

— co <^X<^ co) 

-< 


(13) 


(1 — MJ) j — * 3 ) e ~ v sin a:+- 

Yg (1+2 y) k 2 e~ 2 " sin 2x + j^— (5 + 4y) sin * + 

(4 + 12y + 9y 2 ) sin 3xJ £ 3 e“ w j + ■ ■ • (17) 


where k, the transonic similarity parameter, is 

(t+ l)e 


(14) 


Thus, it is seen that/ is a function of x, y, and the transonic 
similarity parameter k only. In the limiting case, — > 1 

and €— >0, k retains its meaning as a transonic similarity 
parameter. 

Note that in the Prandtl-Busemann iteration procedure 
the order to which the boundary condition is applied at the 
surface of the profile is the same as the order to which the 
iteration has been taken. In the transonic case, however, 
the surface boundary condition is always applied at the 
axis y= 0. In particular, this statement of the boundary 
condition for transonic flow is a rigorous one in the limit as 
both the thickness coefficient and 1— MJ simultaneously 
approach zero but with the same distribution of slope that 
belongs to the family of profiles being treated. A rather 
complete discussion of this point is contained in reference 3. 

Before application of the foregoing considerations to equa- 
tion (10), it should be noted that <j) in that equation is re- 
ferred to the undisturbed stream velocity U. In order to 
introduce c* as the reference velocity, both 4> and the right- 
hand side of equation (10) must be multiplied by U/c*. 

c* 2 

Now, from equation (2), with c=c* and u 2 -\-v 2 =^ji> it 
follows that 


where 


Til M 2 

c y r -*- co 

M* 2 = 

1+1^1 MJ 



(15) 


Comparison of this equation with equation (11) shows that 
the result obtained by means of the Prandtl-Busemann 
iteration method contains implicitly the form of the solution 
required by transonic-flow theory. Moreover, the expres- 
sion for the function f(x,y), namely, 

f(x,y)=—x+(k+-^ g & 3 ) e~ v sin x + 

Yg (1+2 y)k 2 e~ 2v sin 2x+j^— — (5 + 4y) sin x + 

(4+ 12y + 9y 2 ) sin 3xJ fc 3 e“ 3l, + . . . (18) 

satisfies the nonlinear differential equation (12) and the 
boundary conditions, equation (13), to the order k 3 . This 
fact suggests a solution for the flow over the wavy wall, in 
the neighborhood of Mach number unity, obtained directly 
from equation (12). 

SOLUTION OF EQUATION FOR TRANSONIC FLOW PAST 
A WAVY WALL 

Equation (18) suggests the following form for a solution 
of the nonlinear differential equation (12) for transonic flow 
subject to the boundary conditions stated in equation (13): 

f(x,y)= — r + S/» sin nx (19) 

71 = 1 

where the/n's are functions of y only. This form for f(x, y) 
is substituted into equation (12) and repeated use made of 
the following identity: 




co n 
n—0 m = 0 
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When the coefficients of the separate harmonic terms sin nx 
are placed equal to zero, the following system of nonlinear 
ordinary differential equations for /„ results: 


** m— 0 


(m + l)(m+n+l)j m 4 - l/ m + n 4- 1 

7rt=0 

(71=1,2, . . . °°) 


( 20 ) 


Note that equation (20), with the right-hand side placed 
equal to zero, has a solution of the form e~ nv . An iteration 
procedure is then set up in such a manner that the highest 
power of k and of er v appearing in the functions /„ is equal to 
the number designating the iteration step. Thus, the fol- 
lowing sets of iteration equations result: 

Step 1, 

/."-/i=0 


Step 2, 
Step 3, 

Step 4, 
Step 5, 
Step 6, 


//'-4/ 2 =-±/, 2 




/ 3 ,, -9/ 3 =-3/ ] yj 


j " — 4/ 2 = — ^ /1 



3/2/3 

fa ” — 9/3= — 3/1/2 — 6/j/ 4 
/s" — 25/5= - 1 Off,- 1 5/2/3 


/ 2 " - 4 / 2 = — i //- 3 /,/,- 8/2/4 

//' - 1 6/4- - 6/1/3 - 4 /s 2 — 1 0/,/. 

/."- 36 / e = - 1 5 /,/ 5 -|- /s 2 — 24/2/4 


( 21 ) 

( 22 ) 

(23) 

(24) 

( 25 ) 


( 26 ) 


and so forth. 

The right-hand sides of these equations are known functions 
constructed from previously determined quantities in accord 
with the iteration procedure adopted. 

These equations are of the second order, nonhomogeneous 
type with constant coefficients and are readily integrated. 
The resulting expressions for/i, / 2 , / 3 , / 4 , / 5 , and/ 6 , with the 
boundary conditions (13) taken into account, are as follows: 


/,=(l 


11 


256 


k 2 


1861 


3X 256 : 


& 4 ) ke 




2765 



b 33 pY 

! k z e 3v — ^ 

\64^ 

r 64X256 



|_256 1 9X256 2 
65 , 7 


k 2 + 

y+ 


9 


1 


32X256 


8X256 


576X256 ' 72X256 1 
y 3 ^ & 5 e~ w -f- . . . (27) 


Hjn 


16 1 72X256 
4085 


A: 2 


y 


24X256 2 
5 , 815 


432X256 2 
fc2 e -2y_[ ~_ 


/t 4 + 


(r 


11 


1024 

45895 


k 2 


! ) »+( 


L36X256 1 108X256 2 
1 , 11 


P+ 
jfc 2 + 


256 1 288X512 

12245 6245 


144X256 2 72X256 2 y 


128 1 8192 
79 


k 2 

2 y 2 + 


V] 


fcV 


47 


1 


192X256 


12X256 


r 1 23603 /2 

/ 3— 72 "’9.7 VOKR* "■ 


72 1 27X256 
33 


(* 


3X256 
y 4 ^ k s e~ 6vJ r . . . 

259 


(28) 


f 1 | ^ 1.2 

\32 ' 8192 


24 1 72X256 
1765 


) 


fc 2 ) y + 
155 


2 / + 


117 


32X256 


) y2 ] ¥ ' e * V (3X256 2 1 32X256 

+2^6 y 3 ) • • • 0°' 


r 7 4 

390547 ,, . / 

' 7 

32947 A 

_1 536~ 

r 720X256 2 ' 

^384 n 

r 36X256 2 ) 


( 


617 


k 


) y2+ (i( 


11 


128 1 288X256 

/ 35251 2465 

\90X256 2 ‘ 576X256 V 

119 . 25 


96 1 24X256 
455 


2 )y 3 ]i 


48X256 


y 


U (' 


48X256 
25 


96X256 
y 4 ^ Fe- 6 "-f- 


5 2/ 2 + 


(30) 


3840 1 768 
25 


y- 


1536 


y 


384 


y 3 + 


24X256 


y 4 ^ k 5 e 5v T 


(31) 


/s— 


=(- 


91 


91 


23 


13 


V54X2048 1 72X256 
15 


2048 


2048 J 
yi+ 5l20 y ‘) k6e ^ + • 


1024 


y 3 + 


(32) 


Note that the functions _/), / 2 , and / 3 include the terms of 
equation (18), obtained from equation (10) by allowing the 
Mach number to approach unity. 

Equations (27) to (32) may be considered to be essentially 
the nonlinear solution for the flow past a wavy wall of small 
amplitude for stream Mach numbers in the neighborhood of 
unity in the form of a power series in the transonic similarity 
parameter k. Moreover, this solution is identical with the 
one obtained by means of the Prandtl-Busemann iteration 
equations when Mach number unity is approached. 

CALCULATION OF LOCAL MACH NUMBER, CRITICAL 
SIMILARITY PARAMETER, AND 
PRESSURE COEFFICIENT 

GENERAL FORMULAS 

From equation (11), when all terms containing 1 — 
to a power higher than the first are neglected, the expression 
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for the fluid velocity referred to the critical speed of sound 
c* is 

< 33 > 

Now, from equation (2), the relation between q and the 
local Mach number M is 

<Z 2 = — rzii — (34) 

1 +^M 2 

In the transonic approximation, the difference of any Mach 
number in the field of flow from unity is considered a small 
quantity. If terms of only the first power of 1 — M 2 are 
retained, equation (34) yields 

S 2 =l+^fj(M°-l) + . • • (35) 

Hence, from equation (33), 

1-M*=-(1-M.*)!£ (36) 

This equation provides a means of calculating the critical 
value of the transonic similarity parameter; that is, the 
value of the parameter k for which M= 1 at the point of 
maximum velocity on the boundary. 

CALCULATION OF THE CRITICAL VALUE OF k 

For the family of wavy walls of small amplitude (including 
the limiting case of vanishingly small amplitude) with M= 1 
at x—0 and y= 0, equation (36) yields the following relation 
for the determination of the critical value of k: 


(g) -0 


By means of equation (19) together with the expression 
for J ! to / 6 given by equations (27) to (32), equation (37) 
yields the following power series, exact to seven terms, for 
the determination of the critical value of k for the family 
of wavy profiles: 

, . 1 25 73 , 337 74 , 4043 , s , 

^ + 8^ + 384^ + 92l6* + 256X576^’ + 


359381 

270X256 


- 2 F+. . .= 1 


The procedure adopted in order to estimate the critical 
value of k is as follows: From equation (38) the value of 

k can be found for 2, 3, 4, 5, 6, and 7 terms. These values 
of k are, respectively, 1, 0.8990, 0.8644, 0.8504, 0.8424, and 
0.8377. The last two values indicate the approach to the 
asymptotic value of k, that is, the value of k when the number 
of terms in equation (38) is infinite. If the values of k 
approach smoothly to the asymptote, the estimated critical 
value £=0.8377 is very nearly correct. 


Suppose now that both sides of equation (36) are divided 
by [(T+l)e] 2/3 - Then, since 

t (v + 1) 6 


equation (36) can be written as 


1-M 2 

E(tT" l)e ] 2/3 


=— £ _2/3 ' 


The right-hand side of equation (39) is a function of x, y, 
and the parameter k only and is characteristic of the entire 
family of boundary profiles. For the family of wavy walls 
with the critical value assigned to k, equation (39) evaluated 
at the wall becomes 


1-M 2 


KtT l)e] 


273=1.1253(1—0.8536 cos x — 0.0989 cos 2x — 

0.0299 cos 3x— 0.0122 cos 4x — 0.0038 cos 5x — 


0.0017 cos 6r) 


for values of x be- 


Tablc I lists the values of for values of x be- 

[(7+1)4 2/3 

tween — ir and 7 r and figure 2 shows the corresponding 
curve. 

CALCULATION OF THE PRESSURE COEFFICIENT 

Bernoulli’s theorem for a compressible fluid assumes the 
following differential form along a streamline: 


d(p—pj + 2 (>U 2 d 




where 

p pressure in fluid 
p density of fluid 

q speed of fluid 

and the subscript co denotes the quantity in the undisturbed 
fluid. Now, from equation (2) the relation between the 
nondimensional speed q/U and the local Mach number M is 
given by 

*-_ 1+X ill (42) 

772— ! M 2 v**) 

1 + T1 M2 

Ji 

In accordance with the assumption of transonic-flow theory 
that all Mach numbers in the flow differ only slightly from 
unity, equation (42) becomes 

+ . .. (43) 

where powers of 1 —MJ- and 1— M 2 higher than the first 
have been neglected. Then with the pressure coefficient 
defined as 

Cv . m =2=P^ 



] 
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x, deg 


1 — M 2 (y+1) 1/3 

Fioure 2. — Distribution of and - C„, Mcr at the surface for the family of wave-shaped walls. 


equation (41) takes the following form for transonic flow: 


dC„ = d [1 — M 2 — (1 —MJ)\ 
7+1 


After integration, 


C p ,M m = -j-7 [1-M*-(1-M„*)] 

or, with the aid of equation (36), 

^—(1 +|£) (45) 

The right-hand side of this equation can be considered to be 
the first term in a power-series development of C Pt M a in 
1 —MJ. In particular, when the local Mach number first 

attains unity, then ^=0 at x—0, y = 0 and equation (45) 

becomes 


C v , Mcr = ~j(l -M„*) 


a result valid in the transonic range only. Again, if the 
thickness coefficient approaches zero as and M= 1 

at x=Q, y=0, then equation (46) shows that the slope 


dC p , m„ 
d(l-MJ) 
2 


at the critical value M a =M cr ^l is a constant 


— independent of the particular family of profiles 

7+1 

treated. This result is valid whether the approach to Mach 
number unity is made from the subsonic or supersonic region. 

( 7 _|_ 1)1/3 

If both sides of equation (45) are multiplied by — > 


(7 + I) 1 ' 3 „ 

O/Q Tf .M 


-2 k~ m 


0+10 


where the right-hand side depends only on x, y, and the 
parameter k and is characteristic of the family of boundary 
profiles treated. For the family of wavy walls with the 
critical value of k chosen, equation (47) takes the following 
form at the solid surface (?/=0): 

(,y_U 1\U3 

' 2/3 ; C P ,M cr = -2.2507(0.8536 cos z + 0.0989 cos 2x + 
0.0299 cos 32: + 0.0122 cos 4z + 

0.0038 cos 5z + 0. 0017 cos 6x) (48) 

( T _L 1)1/3 

Table I lists the values of - — — C p ,m ct for values of x 
between — ir and w and figure 2 shows the corresponding curve. 
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TABLE I. VALUES OF AND C„ M „ AT 

THE SURFACE FOR THE FAMILY OF WAVE-SHAPED WALLS 


(deg) 

1 — M 2 

(r+lfi'V 

[(t+1)«F 3 

€ 2/3 

-180 

1. 996929 

1. 743190 

-175 

1. 994505 

1. 738342 

-170 

1. 987188 

1. 723708 

—160 

1. 957430 

1. 664191 

-150 

1. 906695 

1. 562722 

-140 

1. 834890 

1.419111 

-130 

1. 743433 

1. 236197 

-120 

1. 634666 

1. 018664 

-110 

1. 510815 

. 770962 

-100 

1. 373602 

. 496536 

-95 

1. 300524 

. 350379 

-90 

1. 224843 

. 199018 

-85 

1. 146916 

. 043163 

-80 

1. 067146 

-. 116375 

-70 

. 903656 

-. 443357 

-60 

. 737110 

-. 776449 

-50 

. 569667 

-1.111335 

-40 

. 404761 

-1. 441145 

-30 

. 250223 

-1. 750223 

-20 

. 119935 

-2. 010799 

-10 

. 031487 

-2. 187694 

-5 

. 007977 

-2. 234715 

0 

. 000007 

-2. 250654 

5 

. 007977 

-2. 234715 

10 

. 031487 

-2. 187694 

20 

.119935 

-2. 010799 

30 

. 250223 

-1. 750223 

40 

. 404761 

-1. 441145 

50 

. 569667 

-1.111335 

60 

. 737110 

-. 776449 

70 

. 903656 

-. 443357 

80 

1.067146 

-. 116375 

85 

1. 146916 

. 043163 

90 

1.224843 

. 199018 

95 

1. 300524 

. 350379 

100 

1.373602 

. 496536 i 

110 

1. 510815 

. 770962 1 

120 

1. 634666 

1.018664 1 

130 

1. 743433 

1. 236197 

140 

1. 834890 

1.419111 

150 

1.906695 

1.526722 

160 

1.957430 

1.664191 

170 

1.987188 

1. 723708 

180 

1.996929 

1. 743190 ; 


n=l, the upper limit r of p goes from 1 to <» and, if n?±l, r 
goes from Oto ®. 

In order to find the recursion formulas for the coefficients 
A"f, the expression for /„ given by equation (49) is sub- 
stituted into the system of differential equations (20). The 
calculation is facilitated by the introduction of the following 
notations: 

71—1 

2 p+n ~2 co 

A;= X vz:A”k* (52) 

3—0 T—v 

and 

(53) 

ff =0 

where the quantities A m ; n arise from the multiplication of 

CO 00 

the two infinite series X A^e~ 2pv and X A%e~ ivv . Note 

p = 0 p = 0 

that the quantities A m ± n are symmetric with respect to the 
upper labels m and n. Then 

/. = **«—*:£ 4; «-*” (54) 

p = o 

When this expression for f n is substituted into equation (20) 
and repeated use is made of the identity 

co co oo n 

SX^.^XS^, 

m=0n=m n=0 m=0 

the exponential terms in y can be eliminated, and the follow- 
ing recursion formulas result: 


GENERAL ANALYSIS 

An examination of the expressions for /, to /„ given by 
equations (27) to (32) shows that the general form of /„ is 

77—1 

2 p+n — 2 co 

fn = j2 e- i2p+n)v x j /' X ^:; i ,+! ' 

V — o 3=0 ? — p 


1 n— 1 

— 2n{Al)' + {AZ)"=— y n X m(n — m)A m '\- m 

4 771=1 

(p=0; n= 2,3, ...«>) (55) 

and 

4 p(p + n)A; — 2 (2 p + n) (A n p ) ' + (A^) " = 


(n= 1, 2, . . . co) (49) 

where, if p = 0, the upper limit of q is n— 1, and, if p ^0, the 
upper limit of q is 2p -j- n — 2. The four-labeled coefficients 
A\ are real numbers calculated from recursion formulas 
obtained from the system of differential equations (20) and 
the boundary condition at the surface of the wavy wall. 
The boundary condition at y= °° is automatically satisfied 
by the form of/„; whereas the boundary condition at the wall 
takes the form 

(//),=«= -S t^l]} (5°) 

Inserting the expression for f n given by equation (49) into 
equation (50) yields immediately the following results: 


and 


Aig=l 


X ( 2 p+n)^S?= X At* 

p = 0 2>=0,1 


(n=l,2, . . . oo) (51) 


where, if n— 1 , the lower limit of p on the right-hand side is 
unity and, if 1, the lower limit of p is zero. Also, if 


1 i i p 

- -ru X m(n—m)A m - n v - m — - c rn X m ( n I m)k 2m A m v lX™ 

4 777=1 £ 777=1 

(z>= 1,2, . . . oo; n=l,2, . . . oo) (56) 

Note that these recursion formulas still contain powers of y 
and k. Both y and k can be eliminated for given values of p 
and recursion formulas containing only the coefficients 
A* f thus obtained. For example, consider equation (55) for 
which p— 0. By repeated use of the identity 


X X ^f= 2 XS.j 

a = 0 (3 = 0 (3=0 a=0 


the following recursion formula is obtained: 


-2n(q+l)A'l 1 0 r +(q + l)(q + 2)8\*A Q % 0 r = 

— y»X X (m + l)(n— m— 1 ) X A^+ t 1 °A n ~r ~ 1 r ~, 

i = 0 m=g—t a=0 


(jl — 2,3, . . . 


°° ; 2 = 0, 1,2, . . . (ti— 2); 


7-=0,l,2, ... no) 


(57) 



8 

where 
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and 


Cn -2 0 

0 a — i 


(2=71—2) 

(2^77—2) 


A number of interesting relations can be easily obtained 
from this recursion formula. Thus, for 2=77—2 and r=0, 

2(71-1) A." x »=y S (7i-m-l)(77i + l)AS^=S 8 A” 1 * 1 8 

"A m— 0 

Then 


A 2 0 1 

Xi -1 o g 

A$ 

a 3 o _jL 

^ 2 0 32 

AI 

a* ° — JL 

^3 0 gg 

• 


25 


48 X 128 


9 

= 5120 


( 58 ) 


From these numerical values the general form is found to be 

A n o _ZL 

“■ 10 7i!4” — 1 

Similarly, for 2=77—2 and r=l, 

n n ~ 2 

A n 0 ^ z A1 o 

^" _1! ( 77 — l)! 4 n - x 1 


In a corresponding manner, more complicated expressions 
can be obtained for •A„t 1 8 when r— 2, 3, . . . . 

Note that in the expression for the disturbance-potential 
function f(x,y), infinite series of the type 


Si/ B 

71=1 


lg-nvfcn + 2 g J n nx A^Li 0 


occur. The ratio of the (77+l)st and 77th terms is 

sin ( 77 + 1)3 A"* 1 0 


Jcye~ 


sm nx 


A n 


_1 

0 

-1 1 


In order for the series to converge, the limit of this ratio as 
77 — > oo must be less than unity. Thus, if the maximum 

value of 7/e~ B (=i) xs inserted and x=0, Cauchy’s ratio 

test yields 

h „ 

(59) 


— 4-<l or 4 
e 4 


Other infinite series occur in the expression for f(x,y) which 
diverge for values of the transonic similarity parmeter k 
considerably less than 4. Thus, consider the recursion 
equation (57) and take 2=0 and r=0. Then 

—2nAi 8+25"y 2 A? 8= 

— (m + l)(n—m — l)A m i 1 ZA n -Z- 1 % 

4 771=0 


A 2 0 

•ft-1 0 — 


8 


C 13 D 943010420 

OXli 0 — &J-12 ol O4l 0 0 

4A18=Ai8+3A8 8+2(A„ 2 8) 2 

10 A? 8=2 Ai 8+10A8 8+15A8 8 A% 8 


27, 


19460 9 4«0J.1lr ISOIOi At 0 l ' (A 3 0\2 

1441i 0 OT 0 1 ^*-do 0 tl 0 oT 2 f-“0 “I 


Also from equation (57), 

4 A! 8= A? 8 

8 A* 8=3^3 8+4 A? 8 A 2 8+3A! 8 
20 Ai 8=6A| 8 + 15A? 8 A8 8+15A? 8 Ag 8 + 10A1 8 

O A6 0 9 46 01 O J20 44 0 ) Q 43 0 43 0 1 O 44 0 42 0_(_C 450 

OA ±2 0 " Al 3 0 I OJTXi 0 -^0 OT & Sll 0 ^0 OT OTli 0 TT 0 01 w«l 0 


and 


30^4| o°= 12^! g+ 15 AS 0 ° A 2 0 0 °+ 10 At 0 °+ 1 5A 2 * A\ o° 

1 O 46 0. 4 46 0 1 Q 43 0 A3 0 1 Q A4 0 42 0 IK 45 0 I 

1 A AI 3 q 4: 0 l y -Tl2 0 -^0 0 A O ^-2 0 ^0 OT ^ Tl 2 0 “T 

q 

O 42 0 A4 0 1 ° C A 3 0\2 
o^li 0 0^2 v^-i oj 

4842 8=20^41 g + 24 Ai g A 2 0 0 0 +l5A* g+274g g AI § + 
24 Aii AH 


The following relations are obtained from the supplementary 
equations (51) : 


O 42 0 A 2 0 

4 TT-q o AA-y 0 

O A 3 0 A 3 0 

o AAq 0 ^1 0 

4A8 8=A18 


K A 5 o 45 o 

O TTo 0 -^l 0 

6A8 8=A?8 


From these relations, together with the ones listed in equa- 
tion (58), 


A 2 0 

-^1 0 

A3 0 

•^*•1 o — 

A 4 0 . 

*“1 0 — 


8 

J_ 

24 

7 

384 


45 0 

+ *-1 0 — 

A 6 0 

o — 


768 

7X13 


72X256 
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and lead to the general rule 

(71=2, 3, ... co) (60) 

From equation (51), tiAs S=A? S; therefore, 

(n=2 ’ 3 ’ ' ' • 00 ) < 61 ) 

where, by definition, 

{3»-5} = lX4X7X10X13X . • - X(3n-5) 

In the expression for the local Mach number in the field of 
flow given by equation (36), the expression bj fjbx occurs. 
This expression, with x=Q, contains infinite series of the 


type 

F=£ nk n e~ nv Ao 2 

(62) 

and 

71 = 1 

G—£ nk n e~ nv A’l □ 

71=2 

(63) 


If the maximum value of e nv {— 1) and the expressions given 
by equations (60) and (61) are inserted, the Cauchy ratio 

4 

test shows that both series converge for &<— • In particular, 

O 

the series expression for F evaluated at the surface (y= 0) 
can bo expressed in closed form. Thus, 


F=k+t : 

71=2 


{371-5} 

ti!4” _1 


k n = 2 — 2 


0 



(64) 


The graph of F against k is a semicubical parabola with the 
4 

cusp point at &=— and F— 2. With the restriction that the 

O 

transonic similarity parameter k be positive and that one 
and only one value of k correspond to a given value of F, 
the permissible values of k and F are confined to the part 
of the parabola lying between the origin (0,0) and the apex 
(4/3,2). The power-series expression for G evaluated at 
the surface can also be expressed in closed form; namely, 


si { 371 5 } 

& “£i(»-l)!4»- 1 



k 


(65) 


This expression, together with the one for F, shows clearly 
that the parameter k cannot be equal to but must be less than 
four-thirds. 

A close examination of the recursion formulas (55) and 
(56) discloses the important fact that each one of the mani- 
fold of power series in k that appear in the functions /„ can 


be expressed in terms of the members of a single dominant 
set of power series. This dominant set consists of one power 
series in k from each namely, the one multiplied by e~ nv 
only. According to equation (49), the members of the 
dominant set of power series are given by 

S n =£ A2° r k n +» (71=1, 2, . . . <*>) (66) 

r= 0 

Several examples are now given to illustrate this impor- 
tant observation. Consider the series that belong to the set 

£ At? k n+2r (71=2, 3, . . . 0°) 

r=0 

From the recursion formulas, 

£ ? k”+*=\ S x 2 

r=0 o 

± At ! jlj S,<+± S.’S.+| S,S.+| S,‘ 


Consider now the set 


Then 


£ Al° r k n+2r (ti= 3, 4, . . . co) 

r=U 

£ A! r ° Si 3 

r= 0 

Zg AX S t *+± St*S a 


Consider the set 

£ Al \ fc n+2r 


(71=1, 2, . . . co) 


From the recursion formula (56), 

£^*" + '=-db s '‘-f s ' s > 

S AS 1 '= ■ -iri-ks s '~h s ' s ' 

S A i ' i ” + *= -Hi s.'-jgj SfS,-^ SfS,- 




256 
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Consider the set 
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X) A J \ k n+2r 


(n=l, 2, ... co ) 


Then 




15 


57 


^AUk 2r+3 2048 512^ 135,2 32 5l25,3 16 




Other examples are 


1 


2048 5l 


S 1 5 - 


S, 3 -S* 


V >41 2 L2r + 1 — 

f=i Zr 8192 ~ 1 256 

S ^ ; *' + '-8U2 s ' s - W W-sW-H s '' s = 


S ^ r 2 * 
r=2 


1 2 7* 2r -f 1 _ 


61 


256X576 


S, s +jI i S,>S',-AsA’+ 


1 0 3 0 1 O O 

256 SlSi 8 5j5s 


The recognition of the existence of a dominant set of power 
series in k represents a major reduction in the complexity of 
the present problem. Thus, consider the array of infinite 
power series contained in equation (66): 


Si=A' 0 o° k+Al ° k*+Al ° 2 ¥+A' 0 ° 3 ¥+ ... ' 
S t =Ai g k 2 +A 2 ° k*+Ao g k°+A 2 0 °k*+ . . . 
S 3 -A 3 8 k 3 +Ai ? k’AAl g k‘+A* 0 ° 3 k*+ ... 
S\=At 0 ° ¥+A* « k*+At g F+ A% ° 3 k 10 + . . . 
S b =A 5 0 8 ¥+A 5 0 ? ¥A-Al ° 2 Jfe»-M8 8 & n + . . : 


An examination of equations (27) to (32) shows that the 
coefficients A% ° of the series S„ appear to be positive and 
monotonically decreasing. The series formed from the first 
column on the right-hand side of equation (67) therefore 
dominates the series formed from succeeding columns. 

co 

Moreover, the first-column series Xj a o 8 k. n has a radius of 

n= 1 

4 

convergence k=— (see equation (61)). The radii of con- 

O 

vergence of the other columnar series therefore are either 
equal to or greater than four-thirds. Similarly, an examina- 
tion of the series S„ shows that the coefficients in each column 
on the right-hand side of equation (67) also seem to form 
positive and monotonically decreasing sequences. This be- 
havior means that Si is the dominant series of the set S n and, 


in fact, of the aggregate of power series in k in the expression 
for the disturbance potential f(x, y) . 

Consider now the series consisting of the first terms of the 
odd-labeled series S x , S 3 , S 6 , ... , that is, 


J2jA 2 \ +10 0 k 2n+1 

71 = 0 


( 68 ) 


According to the theory of power series (and it can be easily 
verified), the radius of convergence of this power series is 
4 

still &=— • Now, a comparison of corresponding terms with 

O 

the dominant series Si shows that 


and 


A 


1 0 
0 1 


11 

256 


>A%° 0 =± 


A 1 0 
■^0 2 


1861 7 

3X256 2 ^ 0 0 15X256 


Thus, if in general A J ° >X 2 " 0 +1 8, then the radius of conver- 
gence of the dominant series Si can be less than the radius 
of convergence of the comparison series given by expression 
(68) and therefore may conceivably be equal to the critical 
value k cr — 0.8377. Moreover, it would then follow that the 
original Prandtl-Busemann small-perturbation method is 
valid for purely subsonic flows only. This conclusion would 
not invalidate other approaches to the transonic-flow problem 
(reference 4). 

Unfortunately, the coefficients Al 8 do not conform to any 
apparent or superficial law, but perhaps a careful study of 
the recursion formulas (55) and (56) and the supplementary 
or boundary relations (51) will yield a rigorous proof of the 
foregoing statements. Otherwise, it remains to calculate a 
reasonable number of the coefficients Al 8 • For this pur- 
pose the development of complete recursion formulas similar 
to equation (57) for the required values of p is worth while. 
Thus, for p — 1, the recursion formula is 

4(n+lM; r ; 1 -2(n+2) (g+ 1 )S A'h 


+(2 + l)(2+2)«"-. 1 'M,? ar i l = 

1 q n-1 -tr 

S (m + l)(n — rn — 

4 t = 0 m—g—t L 


l) 


" A " T _1 , -l + 1+ wi(n — m) S . + 1 A " ' 

s= 0 «=0 


1 




-n(n+l) XMJ M"! 1 r—s 

A 5=0 


(n=l,2, ... co ; q= 0,1,2, . . . n; r= 0,1,2, . . . co) (69) 
where 

6 re — 1 , n 0 

1 —1 


(q=n— 1 or n) 
(q?*n— 1 or n) 


Finally, it may be of interest to give the general formula 
for the Mach number distribution at the surface of the wavy 
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wall. Thus, by means of equations (36) and (49), 
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][/[ 2 — 1 " [ 2 ] m 

l + 2 (n—2m) cos (n— 2m)x A"~ 0 2m £ 

1 1VL CO n= 1 771=0 p= 0 

r (70) 

where — denotes the integral part of n/2. 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., April 26, 1951. 


1. Kaplan, Carl: On Similarity Rules for Transonic Flows. NACA 

Rep. 894, 1948. (Supersedes NACA TN 1527.) 

2. Kaplan, Carl: On the Particular Integrals of the Prandtl-Busemann 

Iteration Equations for the Flow of a Compressible Fluid. NACA 
TN 2159, 1950. 

3. Perl, W., and Klein, Milton M. : Theoretical Investigation and 

Application of Transonic Similarity Law for Two-Dimensional 
Flow. NACA TN 2191, 1950. 

4. Von Kdrm&n, Theodore, and Fabri, Jean: Ecoulemcnt transsonique 

a deux dimensions le long d’une paroiondulde. Comptes Rendus, 
t. 231, no. 23, Dec. 6, 1950, pp. 1271-1274. 


U. S. GOVERNMENT PRINTING OFFICE: 1953 


